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' The Hamiltonian describing a system of strongly correlated electrons coupled to dispersionless 

04 \ phonons was solved numerically for a ring of 8 atoms using the density matrix renormalization group 

■ (DMRG) method. It was found that electron correlation and electron-phonon coupling compete 

against each other, and strong electron correlations suppress the charge-ordered insulating state. 
This allows extended polarons to form in the strong electron-phonon coupling regime. It is shown 
that in this regime two polarons may form a pair via phonon mediated charge-fluctuation interaction. 
Based on the results we propose a novel mechanism of hole pairing. This mechanism could be relevant 
to High-Tc cuprates. 
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I. INTRODUCTION 



i 

. , The role of lattice in the mechanism of the High- Tc superconductivity in cuprates remains controversial. Even 

' though the majority of researchers belieia&.in magnetic mechanisms because of strong-slectron correlations and the 

^ lack of evidence of a strong isotope effectEra, the growing volume of experimental datalj'Q demonstrating the existence 

Q of many anomalous phonon features in cuprates support strong phonon involvement, ft is possible that hole pairs in 

O . cuprates are formed through a cooperative action of strong electron correlations and electron-phonon coupling. The 
purpose of this work is to demonstrate such cooperative effects by numerically studying a model system of strongly 

^ correlated electrons coupled to dynamic phonons. 

^ . In the simplest form such a system is described by the one-band Holstein-Hubbard Hamiltonian: 
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\ where cj^ and b\ stand for electron wiii spin a and phonon creation operators correspondingly, and rii — + ni^. 
0^ ■ Since except for a two-site problemO no exact solution is known for this model, it has been studied mainly by 
numerical methods, but the results obtained so far are not directly apalicable to High-Tc cuprates. The one-polaron 
problem solvpd either by exact diagonalization (ED) on small clustergj'Q, Monte Cario, methodcllj, dynamical mean- 
field theoryllj, or by a special resummation of the strong-coupling perturbation theoryll^l lacks correlation effects. The 
solution of the two-polaron problem obtained by ED for ID and 2D clusters^ also is not relevant to High-Tc cuprates, 
^ ■ since it accounts for correlation effects between just two electrons. To describe these materials realistically one has 
Q ' to include the spin background, i.e. carry out calculations for the system near half-fillins_U» to now such study was 
O ' carried out only in the infinite dimensional limit (dynamical mean-field approximation)ll3~L^. It is not clear if these 
results can be applied to high- Tc cuprates, which are essentially low-dimensional systems. 

Since the amount of calculations grows exponentially with the number of interacting particles, so far this case was 
studied only in approximation of adiabatic {Tllo <^ t) and antiadiabatic (hu S> t) phonons. In the latter -cegime 
^ Hh-h can be transformed by the Lang-Firsov transformation and reformulated in terms of small polaronsE3, but 
■ - ' the results are not applicable to High-Tc cupratesJor which t ^ 1 eV. An interesting modification of the Lang-Firsov 
transformation was employed by Fehske et a/.EZrllj to study the aspects of both adiabatic and antiadiabatic regimes 
simultaneously. 

In the adiabatic regime phonons are treated as "frozen" distortions of the lattice, which can be relaxed into the 
lowest energy state by iterative diagonalizaiion of the electronic part of the Hamiltonian. Such gjjstems studied by 
ED include t-J model (vTO x \ZIO clusters), a two-band version of t-J model {Cuiq0^2 diisterEJ), the three-band 
Hubbard model {CuiO^ clustei — ) and the one-band Hubbard model (VTO x vTO clusteiO). Also few mean- field 
calculations on larger clustersOc3 were carried out. It is clear, though, that the adiabatic approximation misses the 
dynamical aspect of electron-phonon interaction in High-Tc cuprates. 
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The goal of the present work was to study generic effects of electron-phonon coupling in a strongly correlated system 
near half-filling. To this end a small one-dimensional cluster containig electrons interacting with dynamical phonons 
has been studied numerically. We found that a new mechanism of hole pairing based on charge fluctuations mediated 
by phonons arises in this regime. We beleive that it could be a driving force behind the superconductivity in cuprates, 
although a proper study of size and dimensionality effects is required to establish the magnitude of the effect. This is 
left for a further study. 



The Density Matrix Renormalization Group (DMRG) methodEj was applied to a 8-site ring cluster at half-filling, 
with one and two doped holes, and up to 5 phonons per site. The electron-phonon coupling term in ([^) was somewhat 
modified as 



where = (1 — — «ix)' that the term describes coupling of phonons only to a vacant site due to a doped 

hole or charge-transfer fiuctuations of the spin background. Although this definition neglects the difference in the 
site being occupied by one or two electrons, it retains the essential feature of the Frolich term in (|l|): coupling of 
phonons to background charge fluctuations. This satisfies the goal of the study, and at the same time results in a 
sizable decrease of required computations due to lower phonon density. 

The system was studied for a broad range of (7 (4 - 16 eV) and g {0 - 6 eV) parameters, with t = 1 eV and huj fixed 
at 0.1 eV for most of the calculations. This phonon energy was chosen as a characteristic for LO phonons in High- Tc 
cuprates that cause charge transfer frnrn one copper ion to the neighboring copper ion, and show anomalous dependence 
on the momentum and temperaturecZfES. Since in the real system the total number of phonons is not limited, the 
maximum number N^^^^ of 5 phonons per site used in this work is clearly a serious limitation. Unfortunately, the 
increase of N^^^^ results in a rapid growth of the consumed computer memory. The effect of the phonon space 
truncation has been monitored by measuring the average number of phonons per site. For large g/fko values it reaches 
a saturation value, which is ^ 3.3 for N^^^^ = 5 and therefore results for this region are only of a partial merit. On 
the other hand, the most important results (hole pairing regime) was obtained for the small-g sector of the parameter 
space, where effects of phonon truncation are not severe and allow qualitative conclusions be drawn for that part. 



The kinetic energy (t-term in (|I|)) and the total number of holes in the undoped system are shown in Fig. |l| for 
different values of U and g. 



II. METHOD 
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III. RESULTS AND DISCUSSION 



A. Competition between electron correlation and electron-phonon interaction 
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FIG. 1. Kinetic energy and the total number of holes (small symbols) for the undoped system. 
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A rapid decrease in the kinetic energy and a simultaneous increase in the density of holes, defined as a site with no 
electrons {fii = 1), as the value of g is increased suggests that a crossover of the upper Hubbard band (UHB) and the 
lower Hubbard band (LHB) is taking place. As the two bands overlap, the charge transfer excitations create doubly 
occupied sites and the number of holes starts to grow rapidly as illustrated by small symbols in Fig. |l|. For each U 
the overlap g value can be loosely set to the value at which the distinctive bend in the total number of holes in Fig. |l| 
occurs. 

After a sharp decrease the kinetic energy rises again, indicating formation of a state with a smaller degree of charge 
mobility. This happens as UHB sinks below LHB, which corresponds to the charge condensating in the form of 
bipolaron lattice or charge-density wave (CDW) state. This is confirmed by the hole-hole correlation function defined 

by 




shown in Fig. g for the undoped system. 
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FIG. 2. Hole-hole correlation function C(k) for ?7 = 8 eV. k = 1 stands for two holes on neighboring sites, 2 for two holes on 
the next-neighboring sites, etc. Since the total number of holes grows with g, the C(k) is normalized to 1 for each g value 
for illustrative purposes. Vertical lines mark the suggested UHB — LHB overlap region. 

Indeed, as soon as the kinetic energy approaches the minimum, configurations with holes on every other chain 
site (fc = 2 and A; = 4) gain weight, while the weight of the intermediate hole-hole correlation configuration (fc = 3) 
decreases, which reveals the rapid onset of the charge ordering. We like to stress here that the appearance of the 
charge-ordered low-mobility state is delayed up to large values of g in the presence of the on-site correlation energy 
[/, and the delay is proportional to U . In other words, the electron correlations compete against the hole-electron 
interaction and suppress polaron condensation. 

Thus the on-site correlation energy U is effectively reduced by the hole-phonon interaction (|^). Intuitively, the origin 
of this renormalization is easily understood. As a charge fluctuation producing double occupancy on one site and an 
induced hole on another increases the energy of the systemj-by U , it also couples to the phononic field via (H) thus 
reducing the effective U to Uef f < U (for discussion see alsollj) . The minima of the kinetic energy curves correspond 
to the maximum UHB — LHB overlap, and therefore Ueff — 0. The relation between the effective Hubbard energy 
Ueff(U,g) and g depends on the g/Tiuj ratio and is nonlinear, since the increase in g will cause the increase in the 
average number of phonons per site and, consequently, the g dependence. However, as can be seen from Fig. 0, the 
value of U corresponding to Ueff = increases proportionally to g, i.e. the relation is approximately linear for this 
system. We beleive that this is an artifact of the Np^^^ truncation as discussed above, since it has been found that 
the phonon density saturates shortly before the kinetic energy minimum, and the regime at the minimum and above 
does not allow the number of phonons in the system to grow along with g. 
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FIG. 3. Kinetic energy and the total number of 



for the system with 0, 1 and 2 doped holes {U — 8 eV). 



B. Nature of the doped hole: Extended polaron formation 



We now turn to the system doped with one hole. One can notice in Fig. that the kinetic energy of the one-hole 
doped system increases slightly as hole-phonon interaction is turned on, suggesting localization tendency of holes, 
apparently due to phonon drag. The extent of phonon dressing is illustrated by the hole-phonon density correlation 
function P(k) defined by 

where n^^j, stands for the number of phonons at the site i + k. P(k) for the system with one doped hole is shown in 
Fig. ^. One can see that for small g values phonons concentrate around the doped hole and the phonon cloud shrinks 
rapidly as g increases. Since the kinetic energy increases at the same time, it is most likely that a small polaron is 
forming. As the value of g is further increased, approximately half-way to the UHB — LHB overlap, a phonon cloud 
surrounding the doped hole starts to spread again, so that by the time the system reaches the overlap, the cloud is 
extended over the whole cluster. This is accompanied by a decrease in the kinetic energy, so that this result suggests 
that an extended polaron is formed as a precursor for the transition into a high-mobility phase at the UHB — LHB 
crossover. 
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FIG. 4. Hole-phonon density correlation function P(k) for U 



eV. The Pik) is normalized to 1 for each g value. 
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This behavior is very different from that of a conventional Holstein model. It is well known that for the conventional 
model strong electron-phonon coupling {g/Twj > 1) results in a polaron band collapse (self-trapped polaron) for a 
single electron, and the CDW formation at half-filling. Here, on the other hand, strong electron correlation competes 
against the electron-phonon coupling, allowing a phase in which charge becomes mobile before the polaron ordering 
happens - possibly the extended polaron phase. It is interesting that a similar effect of electron correlations enhancing 
the conductivity and causing an IM transition was recently found for a disordered Hubbard model in two dimcnsionsE3. 



We now discuss the system with two doped holes. It will be shown that in the range of parameters that corresponds 
to the extended polaron phase, the two doped holes bind into a pair larger ly due to reduction of the background 
charge fluctuations existing in the phonon cloud surrounding the pair. By defining the hole-pair binding energy as 



where Eg{n) is the ground-state energy for the system doped with n holes, one can study the hole-pair formation 
and estimate the pairing energy. The obtained binding energy is shown in Fig. ^ as a function of g and U . One can 
see that the {U , g} parameter space can be divided into four regions. For small values of g Ef, is positive, i.e. no 
hole-pairing occurs. As g increases, Ei, becomes negative if U is large enough for the gap between UHB and LHB to 
exist, and a hole-pair is formed. As g increases futher to the point where the UHB and LHB overlap, the binding 
energy becomes positive again and the pair formation becomes unfavorable. Finally, as the system crosses over into 
the charge-ordered region, Eb becomes negative again. Therefore, there are two regions of hole-hole attraction — one 
in the charge-ordered state, and another in the vicinity of the UHB — LHB overlap. The first one obviously results 
in aggregation of doped and induced holes into the charge superlattice. The second one corresponds to the extended 
polaron regime where the Hubbard gap just closes and the polarons become unstable to pair formation. Clearly, as 
applied to cuprates the region of interest is the second one, as it has a possibility of resulting in superconductivity. 



C. Charge-fluctuation mechanism of polaron pairing 



Et^Eg{2) + Eg{0)-2Eg{l), 
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FIG. 5. Hole pair binding energy Eb- 
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FIG. 6. Components of the hole pair binding energy for U — S eV. The binding energy is shown by the dashed curve. 

To understand the origin of the binding force in this region we separate Ei, into "polaronic" {t and g terms), 
"phononic" (hu; term), and "correlation" [U term) inputs, by calculating expectation value of the corresponding 
terms in the studied Hamiltonian for the ground state. The result for J7 = 8 is shown in Fig. |6[ Two important 
conclusions follow. First of all, one can see that the correlation energy plays a crucial role in the pair formation 
mechanism, giving most of the "binding" (negative) input into Eh. On the contrary, the effect of dynamic "polaronic" 
interaction is to the most part negative or "anti-binding" . A peculiar dip in the "polaronic" curve and corresponding 
bumps in "phononic" and "correlation" curves at g ~ 0.5 eV correspond to the small polaron — extended polaron 
transition mentioned above. In fact, for larger U values the dip in the polaronic binding energy creates a dominant 
input into the hole-pair binding for small values of g, i.e. before the exteded polaron is formed, which corresponds to 
a conventional small polaron binding mechanism. 

Secondly, the behaviour of Ei, components clearly correlates with the UHB — LHB crossover, when the extended 
polaron becomes unstable. As g increased through the crossover, the polaronic component of Ei, remains positive 
and increases, while the correlation component stays negative and decreases. As the two bands begin to overlap, the 
polaronic and the correlation parts reach their maximum and minimum correspondingly. When the overlap is at its 
largest (C/e// — 0), the two components change signs and diverge as soon as the system reaches the charge-ordered 
regime. 

From the definition of Eh we see that it compares the energy of two separate doped holes and their backgrounds 
with the energy of two holes confined into one background plus the undoped background itself. Therefore, the positive 
sign of the Eb polaronic component, for example, signifies that the system with two separate doped holes allows more 
hole motion than the system with the two holes bound together. By analyzing Fig. ^ along these lines we propose 
the phonon mediated charge-fluctuation mechanism of hole pairing as described below. For illustration purposes the 
model is generalized to two dimensions. 
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FIG. 7. Two regimes of doped holes binding: (a) small polarons, (b) extended polarons, (left) unpaired state, (right) paired 
state. Large open circles stand for doped holes, small black and white circles for background charge-transfer fluctuations, while 
shaded regions represent phonon clouds. 



We start from the state with .g = and U large enough for the Hubbard gap to be open. Since two doped holes 
confined to one region strongly interfere with each other's motion Eb is positive. As the hole-phonon interaction is 
turned on, conventional small polarons are formed, and motion of a hole is reduced due to renormalization of its 
mass. In this regime two holes can combine their phononic clouds, thus increasing their mobility as compared to the 
separated holes (Fig. ^) . This is a regular polaron coupling mechanism, and it applies to the hole-pairing regime for 
large U and small g mentioned above. As the g/huj ratio is increased the phonon cloud of a single polaron spreads on 
neighboring sites forming an extended polaron (Fig. ^) . The role of correlation interaction U becomes crucial here 
since it delays the charge-ordered state formation and allows extended polarons to form. In this regime holes tend 
to share their phonon clouds, which will decrease their kinetic energy relative to the unpaired holes. Altogether this 
explains the dip observed in the polaronic component of Ei,. 

The main reason for extended polarons to share their phonon clouds is the charge-transfer fluctuations of the 
background electrons. Extended phononic field causes charge-transfer background fluctuations (|, I) ^ (0, ti) via 
the hole-phonon interaction (see Fig. |^b). The double ocupancy will cost on-site correlation energy U (in addition to 
the hu! energy of the phononic fleld), and unless the UHB and LHB overlap, the two polarons prefer to share their 
domains to minimize the number of the (0, ti) pairs. This is corroborated by large negative values of the correlation 
(U) component of Ef, in Fig. |^. On the other hand, as the bands overlap, the "fluctuation" pair formation costs less 
and less energy, so that the tendency for polarons to stay together will decrease, until the polaron-pair formation will 
become unfavorable {Ej, becomes positive again in Fig. ra). This phonon- induced corxclation mechanism of pairing is 
closely related to the idea of " correlation bag" mechanism proposed by GoodenoughS. 
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FIG. 8. Ising and charge-transfer components of the correlation component of Eb for U ■ 
by the dashed curve. 
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Since the Hubbard "correlation" component U includes both spin-fluctuation (as in the t — J model) and charge- 
fluctuation parts, it is important to estimate their effect on pair formation separately. To this end we evaluated the 
"t — J" part by calulating the ground state expectation value of the Ising term: ^(X^i ^t^i+i)- The charge part 
was defined as the difference between the U term in Fig. 6 and the spin fluctuation contribution. The resulting 
components are shown in Fig. H. It can be seen that the spin contribution is dominant in the small-polaron regime 
as expected, but it becomes weaker as the UHB - LHB crossover occurs. Instead, the charge-fluctuation input into 
-Ef), which is suppressed in the small-polaron regime, quickly becomes dominairt as the Hubbard gap decreases and 
the charge-transfer process is activated. 
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D. Relevance to the High- Tc cuprates 




g 



FIG. 9. A generic U-g pliase diagram for tlie one-band Hubbard-Holstein model. PSC stands for polaronic superconductivity, 
and CFSC for phonon mediated charge-fluctuation superconductivity. 

A schematic U-g phase diagram in Fig. |^ summarizes the qualitative picture obtained in this work in the form 
of suggested behavior of a real physical system (a high- Tc cuprate, for example). The relative position of the UHB 
and the LHB determines its electronic behavior. When the Hubbard gap is open (large U, small g), the system 
is a doped insulator with polaronic conduction. When U is small, or the effects of U and g are competing, the 
UHB and LHB overlap and the system is metallic. For small U and large g the insulating bipolaronic-lattice/CDW 
regime prevails. The superconducting transition happens as polarons become unstable to pair formation. For large 
U values the pairing mechanism is the conventional polaronic coupling, while for lower U values in the vicinity of 
the polaron-metal crossover the mechanism is that of phonon mediated charge-fluctuations described above. While 
the superconducting transition temperature in the polaronic regime is expected to be low due to the strong mass 
renormalization, in the polaron-metal crossover regime the spread of a polaron would result in a much smaller mass, 
and thus could result in a high transition temperature. 

Qualitatively it is clear that the proposed charge-fluctuation mechanism will produce the maximum pair-binding 
energy at the Insulator-to-Metal (IM) crossover, as seen in Fig. |6[ Indeed, by suppressing the superconductivity Ly 
magnetic field it has been shown that the IM crossover happens at the same doping as the maximum Tc is reached lj. 

The doping driven IM transition in LSCO suggests a percolative origin of conductivity at low temperatures. In fact, 
such a mechanism involving percolation between larse. polarons existing at low temperatures was proposed to explain 
IM transition and CMR phenomenon in manganitesL^lS. The present results demonstrate that the large (extended) 
polaron regime can arise naturally in doped Mott-Hubbard insulators with strong electron-phonon coupling. As the 
number of doped holes increases extended polarons overlap, and the percolative conductance developes. Furthermore 
it was observed for manganites that local dynamic coupling of hole and phonon (vibronic state) develops at the IM 
transitionEj. In cuprates such local dynamic coupling could lead to superconductivity of extended polarons. The 
anomalous behavior of LO phonons recently observed by neutron inelastic scattering could be the manifestation of 



By solving the one-band Hubbard hamiltonian coupled to dynamic phonons by a Frolich-type term near half- 
filling several novel observations have been made: 1) Electron correlation and electron-phonon coupling compete with 
each other. Thus the charge-ordered insulating state is suppressed by electron correlations, and extended polaron 
formation becomes possible for strong electron-phonon coupling. 2) The extended polarons formed by doped holes 
bind into pairs, and their binding energy is strongly enhanced by electron correlation contribution. Based upon these 
observations we propose that phonon mediated pairing interaction based upon charge fluctuations in the half-filled 
system can be a relevant mechanism of superconductivity in High- Tc cuprates. 
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IV. CONCLUSION 
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